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We report diffusion Monte Carlo calculations on the phase diagrams of para-}i2 and ortho-Y)2 
adsorbed on top of a first layer of the same substances on graphene. We found that the ground 
state of the second layer is a triangular incommensurate solid for both isotopes. The densities 
for promotion to a second layer and for the onset of a two-dimensional solid on that second layer 
compare favorably to available experimental data in both cases. 
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I. INTRODUCTION 

Graphite and graphene are closely related forms of car- 
bon in which the atoms are located in the nodes of a hon- 
eycomb lattice. The main difference is that graphite is 
formed by a whole stack of the two-dimensional carbon 
sheets that constitute a single graphene layer J^ Adsorp- 
tion of gases on top of graphene and graphite are ex- 
pected to show similar trends. This is in fact what has 
been shown by computer simulations of the phase dia- 
grams of adsorbed quantum gases on graphene, both in 
the first, '^"^ and second^^ layers. In general, those calcu- 
lations render phase diagrams which are very close to the 
experimental ones on graphiteJ^^— The only appreciable 
difference is the binding energy of the adsorbate species 
on top of the carbon structure, bigger in the graphite 
case. Unfortunately, to our knowledge, there are no ex- 
perimental data yet on adsorption of quantum gases on 
graphene to be compared to. 

Most of that work, both from the experimental and 
theoretical sides, has been devoted to the determination 
of the phase diagram of the first layer of quantum gases 
and solids adsorbed on top of graphite.™ However, there 
has been some measurements of the properties of a sec- 
ond "^He sheet on top of an incommensurate helium layer, 
directly in contact with a graphite surface jiiii^ii^"— work 
that can be compared to the simulations on the same swh- 
jecti ^^d^ On the other hand, the second layer of molec- 
ular hydrogen on top of graphite has been less explored. 
Experimentally, the second layer of molecular hydrogen 
and deuterium adsorbed both on graphite and MgO has 
been studied using calorimetric measurements. ^"^^'^ In 
particular, Ref. [23 studies the promotion to a second 
layer and the phase diagram of pure para-H2 (P-II2) 
and ortho-D2 (0-D2) second layers on graphite. One of 
the main conclusions of that work is that the ground 
state of both isotopes in the second layer is a quasi two- 
dimensional solid. Those calorimetric measurements sug- 
gest triple points at T ^ 6 K for H2 and T ^^ 11 K for 
D2. For lower temperatures and densities those solids 



seem to coexist with infinitely diluted gases. The main 
goal of our present work is to perform diffusion Monte 
Carlo calculations to determine the phase diagram of the 
second layer of H2 adsorbed on a first layer of H2 on 
graphene, and of the second layer of D2 adsorbed on the 
first layer of D2, also on graphene. The results so ob- 
tained will be compared to the experimental ones for the 
same systems on graphite?^ the only ones available, in 
order to assess the possible differences. In particular, we 
try to ascertain if the oblique structure suggested for D2 
in Refs . |2 3l, from neutron diffraction experiments, and 
in Ref. [2J, from low-energy electron-diffraction (LEED) 
measurements, is more stable than an arrangement con- 
sisting of two incommensurate triangular solids of differ- 
ent densities. 

In the next Section, we describe the theoretical method 
used in the microscopic study of the adsorbed phases. 
The results obtained for both P-II2 and 0-D2 on graphene 
are shown in Sec. III. Finally, the main conclusions are 
discussed in Sec. IV. 



II. METHOD 

The diffusion Monte Carlo (DMC) method is a stochas- 
tic technique that allows us to obtain the ground state 
of a zero-temperature system of bosons. Since the 
species adsorbed (j>-H2 and 0-D2) on graphene are both 
bosons, we can obtain through DMC the arrangement of 
molecules with the lowest energy for each surface den- 
sity. In order to reduce the statistical variance of the 
many-body problem, the algorithm uses a guiding wave 
function $ which enhances the occupation probability 
in places where the hydrogen density is expected to be 
large. ^^ In general, $ depends on the coordinates of all 
atoms or molecules in the simulation cell. However, in 
this work, we will consider that the carbon atoms of the 
graphene layer are kept in fixed positions, what means 
that their overall effect on the hydrogen molecules can 
be described as an external potential. Within this ap- 



proximation, the guiding wave function depends only on 
the positions of the N hydrogen molecules (ri,. . .,rjv). 
We chose, 

$(ri,...,rjv) = $j(ri,...,rjv)$i(ri,...,rArJ 

X $2(rjVi+i,---,riv) , (1) 

where Ni {N2 = N — Ni) is the number of molecules 
adsorbed in the first (second) layer, and $j is a Jastrow 
wave function used to take into account the H2-H2 and 
D2-D2 correlations,'*' - In particular. 
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r,j being the distance between the center of mass of hy- 
drogen molecules (considered as spheres and interacting 
through the isotropic Silvera and Goldman potential^). 
This potential was developed to describe the hydrogen- 
hydrogen interaction in bulk and it was successfully used 
to reproduce the phase diagram of the first layer of both 
H2 and D2 on graphenei^ii The variational parameter b 
in Eq. © was taken to be 3.195 A for both H2 and 
D2, in agreement with variational optimizations done in 
previous simulations of similar systems i^ii 

The purpose of the remaining terms in Eq. (dJ, $1 and 
$2, is to describe the localization of hydrogen molecules 
in the first and second layers on top of graphene, re- 
spectively, and also the correlations due to the carbon- 
hydrogen interactions. Concretely, 
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$„(ri,r2,...,rAr„) = 1111^^? 
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JJexp(-a„(zj - z„)^) 



(3) 



where (a;sitc,2/sitc) are the crystallographic positions of 
the solid lattice. The c„ parameters were taken to be 
the same for n= 1,2, i.e., a linear interpolation between 
the values corresponding to densities in the range 0.08 
A-2 {cn = 0.61 A-2) and 0.10 A^^ (c„ = 1.38 A^^) 
for R2r and between 0.08 A^^ (c„ = 1.11 A-^) and 
0.11 A-2 (c„ = 2.93 A-2) in the case of Dj.^ If the 
hydrogen density within the considered layer was not in 
those ranges, we used the linear extrapolated c„ value. 

To model the second layer of molecular hydrogen on 
top of a first layer of the same substance, we followed 
closely the prescription of Ref. M, in which a second layer 
of ^He on graphene was simulated. Basically, for a fixed 
total hydrogen density, we considered only the arrange- 
ment for which the total energy per molecule was lower. 
In practice, this means that we have to take a fixed solid 
density for the first layer, and to change the number of 
molecules (if we have a liquid) or the lattice constant (if 
we have a solid) in the second layer. The structure of the 
second layer solid was assumed to be incommensurate 
with respect to the one on the first layer, i.e., the phase 
diagram was assumed to be of the same type as that of 
helium on graphene. However, to verify that this was so, 
we considered also the oblique commensurate structure 
proposed for D2 from neutron diffraction^^ and LEED^l. 
We checked the stability of that structure for both iso- 
topes. To be able to treat incommensurate second lay- 
ers, we used different periodic boundary conditions for 
the first and second layers. No exchange of molecules 
between the first and second layer was allowed. Impor- 
tantly, we did not fix the positions of the molecules clos- 
est to the graphene surface, i.e., we took into account the 
zero point motion of all the hydrogen molecules. For both 
isotopes, we considered hydrogen densities up to those of 
promotion to a third layer, obtained experimentally for 
the same systems on graphite 1^ 



with n = 1, 2 for the first and second hydrogen layer, 
respectively, and Nn the number of molecules in each 
layer. Vij represents the distance between the center or 
mass of each molecule, i, and each of the Nc carbon 
atoms, J, in the graphene layer. Each hydrogen molecule 
interacts with each of those carbon atoms by a potential 
of Lennard- Jones type, whose parameters are taken from 
Ref. [23. This part of the guiding wave function is again 
similar to the one used in previous works to describe the 
first layer of H2 (Ref. 0) and D2 (Ref. 0) on graphene, to 
the point that we took the same oi and zi parameters as 
in those works. Thus, ai = 3.06 A~^ for H2 and ai = 5.2 
A"2 for D2. be ^ 2.3 Aand 21 = 2.9 Afor both hydrogen 
isotopes. For hydrogen molecules on the second layer, 
be was kept constant, and 02 and Z2 were variationally 
optimized. The optimal values were 02 = 1.53 A^^ for 
both H2 and D2, and ^2 = 6 A(H2) and 5.8 A(D2). 

When the phase to be described is a quasi two- 
dimensional solid, $„ is multiplied by a Nosanow term, 

JJ exp{-c„[(xi - a;sito)^ + ivi - 2/sito)^]} , (4) 



III. RESULTS 



H2 



The phase diagram of the second layer of H2 on 
graphene can be extracted from the data displayed in 
Fig. [1] The full squares correspond to the energy per 
molecule of a triangular incommensurate solid adsorbed 
on a single layer on top of graphene, and were taken from 
Ref. y. When we considered only a hydrogen layer, in- 
commensurate means that there is no registry of the hy- 
drogen molecules with respect to the carbon sheet. The 
third-order polynomial fit displayed on top of them is a 
guide-to-the eye. To study second layer structures, we 
put on top of one of these incommensurate structures a 
set of hydrogen molecules described by a guiding wave 
function with C2 = (|4]). The dimensions of the sim- 
ulation cells were determined by the density of incom- 
mensurate solid in the first layer since, in all cases, the 
simulation cells comprised 120 H2 molecules close to the 



graphene surface. Then, we put enough hydrogen on top 
of them to produce surface areas (the inverse of density) 
in the range displayed in Fig. [T] To be sure that the en- 
ergy per molecule considering both hydrogen sheets was 
the minimum one, we performed several sets of simula- 
tions with different incommensurate solid densities (in 
the first layer). In particular, we considered 0.090, 0.095, 
0.10 and 0.105 A~^. Our results indicate that on the 
second layer, the arrangement with lowest energy per H2 
molecule is a liquid on top of a solid whose density is 
0.095 A~^. Those results are displayed in Fig. [T]as full 
circles. The lowest limit of the surface area corresponds 
approximately to the experimental value for H2 promo- 
tion to a third layer (5.80 A^)^ 

We modeled the second layer incommensurate trian- 
gular solid on the same principles, i.e., we considered the 
same densities as for the liquid case in the first hydrogen 
sheet, and distributed the atoms on the second layer as 
to produce the total densities displayed as surface areas 
in Fig. [TJ Here, incommensurate means that the second 
layer is registered neither with respect to the first layer 
nor with the underlying graphene. The minimum ener- 
gies per H2 molecule corresponded to an arrangement in 
which the lower layer density was 0.10 A~^; those data 
are displayed in Fig. [T]as open squares. The line on top 
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FIG. 1: (Color online) Maxwell construction (straight line) to 
determine the limits of the phase diagram of a second layer of 
H2 on top of the same substrate on graphene. Full squares, 
energy per H2 molecule on a first layer triangular incommen- 
surate solid; full circles, same data for a second layer liquid on 
a triangular solid; open squares, a triangular incommensurate 
solid on top of first layer solid. The open triangle represents 
the energy per H2 molecule of the obfique structure proposed 
in Ref. [231, while the full triangle is the energy of a second 
layer 4/7 structure. The curves are three-order polynomial 
fits to the simulation results. 



of them represents also a third-order polynomial fit. We 
considered also two second layer commensurate (with re- 
spect to the first layer) solids: a 4/7 lattice;^ and an 
oblique bilayer structurCj^- both proposed originally to 
describe second layers of '^He on graphite. The latter 
was suggested to be stable for the second layer of H^^ 
and D2 —1^ on graphite. The energy per H2 molecule of 
those registered phases is represented by a full and open 
triangle, respectively. We can see that both structures 
are metastable with respect to a set of two incommensu- 
rate layers, since their energies per molecule are larger. 
Their corresponding energies per II2 molecule are -274.9 
± 0.1 K (oblique structure), and -304.6 ± 0.1 K (4/7), 
versus -293.0 ± 0.1 K and -309.38 ± 0.08 K of the cor- 
responding incommensurate structures of the same den- 
sities (0.173 and 0.157 A'^, respectively). A 7/12 com- 
mensurate solid, also proposed to be stable for helium,— 
and not displayed for simplicity, was also considered and 
found to be unstable with respect to a second layer in- 
commensurate solid. 

With all of that in mind, we can draw a double-tangent 
Maxwell construction to obtain the phase diagram of the 
second layer of H2 on graphene. The slope of that line is 
minus the internal pressure at which the transition takes 
place, and has to be positive for stable arrangements. In 
addition, if several transitions are possible, one has to 
consider only the corresponding to the lowest pressure. 
The straight line displayed in Fig. [I] that joints a single 
layer incominensurate solid and its second layer counter- 
part fulfills all the necessary requirements. We can see 
that the energy per II2 molecule of a second layer liq- 
uid (full circles) is always above the double-tangent line. 
This means that a second layer liquid is unstable with 
respect to a mixture of a single layer solid and another 
with two incommensurate sheets. The surface areas at 
which the slopes of both three-order polynomial fits are 
the same define the stability limits for the phases involved 
in the transition. In our case, those correspond to an up- 
per density limit for a single layer solid (equivalent to a 
second layer promotion) of 0.094 ±0.002 A^^, and a low- 
est density value of 0.160 ± 0.002 A^^ for a two layered 
solid. Both results are in excellent agreement with the 
calorimetric values of II2 on graphite:^'^ 0.094 and 0.165 
A~^. The chemical potential for II2 at the monolayer- 
bilayer transition derived from our simulations was -161 
± 2 K. This means that a second layer of H2 is still stable 
with respect to the formation of a bulk H2 crystal, whose 
ground state chemical potential is -92.3 K— . This also 
means that the second layer solidifies at H2 densities as 
low as 0.060 A~^. As in "'Hcj^iiiii^ the first layer suffers a 
compression that produces a ~ 6% increase of its density 
upon the adsorption of the second hydrogen sheet. 

We also compared our simulation results for H2 on a 
second layer to similar data obtained for a pure two di- 
mensional (2D) system. In Fig. [2l we show the energy 
per II2 molecule as a function of density, but only for the 
molecules on that second sheet. This energies are taken 
from the same set of simulation results as the ones dis- 



played in Fig. [1] i.e., the molecules of the first layer are 
not kept frozen. A glance to the previous figure indicates 
that the binding energies of these molecules are smaller 
than the ones located in the first layer by at least a factor 
of two. This is true for both liquid (full circles) and tri- 
angular incommensurate (full squares) phases. The lines 
on top of each set of symbols are least-squares fits to the 
expression 



E/N - {E/N)o + a{p - pof + b{p - pof 



(5) 



where p is the hydrogen density in that second layer and 
po stands for the density at which the energy per H2 
molecule has a minimum {{E/N)q). The parameters ob- 
tained for the liquid and solid phases, together with their 
first layer (taken from Ref. 4) and pure two-dimensional 
(from Ref. l29t ) counterparts, are shown in Table HI Those 
results indicate that a pure two-dimensional system of 
H2 molecules is a worse description for the second than 
for the first layer of H2 on graphene, and that the 2D 
equation of state is worse for a liquid phase than for a 
solid one. For instance, the energy differences, once sub- 
tracted the infinite-dilution energy (i?oo), between a 2D 
system and a first layer of H2 on graphene are '^ 6.5 % 
for a liquid, and ~ 5 % for a solid. In comparison, those 
same differences between a 2D and a second layer system 
are ~ 71 and 22 %, respectively. The same, but to a 
lesser extend, can be said of the differences between the 
values oi pq: ^ 2.3 % for the 2D and first layer solids 
difference, gap that increases up to ^ 4.2 % if instead of 
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FIG. 2: (Color online) Energy per H2 for molecules located on 
the second layer. Full circles are the results of a liquid phase, 
and full squares, for a triangular incommensurate solid one. 
The error bars are similar in both cases and for simplicity 
they are only shown in the liquid case. Lines are three-order 
polynomial fits to the respective data. 



TABLE I: Density and energy per molecule for the phases 
displayed in Fig. [2] from fits using Eq. ^. The results for a 
pure two dimensional system (2D, Ref. i29i) and a first layer 
of H2 on top of graphene (gr, Ref |j) are also shown for 
comparison. Eoo is the energy per H2 molecule in the infinite 
dilution limit for each system, {E/N)o is the minimum energy 
per II2 molecule and po is the density that corresponds to that 
energy. 



Liquid 



Solid 



2D 



Soo (K) 

iE/N)o (K) 

iE/N)o - E^ (K) 





-21.43 ± 0.02 

-21.43 ± 0.02 

0.0633 ± 0.0003 



-23.453 ± 0.003 
-23.453 ± 0.003 
0.0673 ± 0.0002 



gr 



E^ (K) 

iE/N)o (K) 

{E/N)o - Soo (K) 



-431.79 ± 0.06 

-451.88 ± 0.03 

-20.09 ± 0.07 

0.05948 ± 0.00005 



-454.1 ± 0.3 

-22.3 ± 0.3 

0.0689 ± 0.0006 



second layer 



E^ (K) 

iE/N)o (K) 

{E/N)o - E^ (K) 

PO (A-^) 



-114.4 ± 0.6 

-126.9 ± 0.5 

-12.5 ± 0.8 

0.055 ± 0.001 



-133.7 ± 0.3 

-19.3 ± 0.7 

0.0646 ± 0.0009 



a first layer we have a second one. Those same percent- 
ages grow to ~ 6.5 % and 15 % for the same comparisons 
for liquid phases. Moreover, Table U gives us another in- 
teresting piece of information: the upper solid is more 
stable with respect to a liquid arrangement than in a fiat 
structure or a first layer sheet. This is because the the 
energy difference at the respective equilibrium densities 
is highest in the second layer case. The corresponding 
values are 6.8 K (second layer), versus 2.023 K (2D), and 
1.8 K (first layer). We can also see that the density min- 
imum for the second-layer solid phase (0.0646 ± 0.0009 
A~^) is comparable to the equilibrium density mentioned 
above (0.060 ± 0.002 A'^), obtained from the Maxweh 
construction including the whole system. 



D2 



This subsection will closely mirror the previous one 
since we studied the second layer of D2 on graphene fol- 
lowing the same steps. Our results are summarized in 
Fig. [3l where all the symbols and lines have a similar 
meaning to those of Fig. [TJ The only difference is that 
the structures with the minimum energy per deuterium 
molecule are those whose first layer density is 0.105 A~^ 
for both the second layer solid and liquid phases. To be 
sure of that, densities in the range 0.095 to 0.110 A~^ 
were tested. The first layer results (full squares) are now 
taken from Ref. 0, and the triangles represent the same 
registered phases suggested above for H2 . As one can see, 
these commensurate phases are still unstable with respect 
to a set of two incommensurate solid deuterium layers. 
As in II2 , we have also that the energies per D2 molecule 



for a second layer liquid are above the double-tangent 
Maxwell construction. This means that the phases in 
equilibrium are again a first layer incommensurate trian- 
gular solid of density 0.100±0.002 A~^ and a second layer 
incommensurate solid whose total density is 0.175±0.002 
A^^. Both results are again in excellent agreement with 
the calorimetric data of Ref. [23 on graphite: 0.099 A^^ 
for second layer promotion, and 0.178 A~^ for the mini- 
mum density at which the double solid is stable. Since in 
this last case, the density of the lower layer is 0.105 A~^, 
we can state that there is also a compression of the first 
layer of around a 5 % upon the adsorption of a second 
layer of D2 on top of D2 . As in the case of H2 , our simu- 
lation results allow us to calculate the chemical potential 
of D2 at the monolayer-bilayer transition. The result was 
-180 ± 2 K, a healthy 30% larger than the corresponding 
to bulk D2. We found also that both the oblique com- 
mensurate structure and the 4/7 one were unstable with 
respect to the double incommensurate arrangement, as 
can be seen in Fig. ([3]). The values of their energies per 
D2 molecule were -326.1 ± 0.1 K and -352.7 ± 0.1 K, 
versus -341.1 ± 0.1 K and -358.0 ± 0.1 K for the double 
incommensurate for the same densities (0.186 and 0.164 
A~^, respectively). 

We can also study the second layer by itself and com- 
pare the results to those of a first layer of D2 adsorbed on 
graphene, and to a pure 2D system. That can be done 
with the help of Fig. |4] and Table |TT1 The conclusions 
we can draw from this set of information are similar to 
those already described in the H2 case: the equilibrium 
densities are essentially compatible with each other, even 



TABLE II: Density and energy per D2 molecule for the phases 
displayed in Fig. [J] obtained by the same means as those of 
II2. The results for a pure two-dimensional system (2D, Ref. 
[291 . with only data for the solid phase) and a first layer of D2 
on top of graphene (gr, Ref. [3 for the liquid and this work of 
the solid) are also shown for comparison. The variables have 
the same meaning as in Table [D 



Liquid 



Solid 



2D 



E^ (K) 

{E/N)o (K) 

{E/N)o - £^00 (K) 



-42.305 ± 0.005 
-42.305 ± 0.005 
0.0785 ± 0.0002 



gr 



£00 (K) 

{E/N)o (K) 

iE/N)o - E^ (K) 

PO (A-^) 



-464.87 ± 0.06 

-497.2 ± 0.9 

-32.3 ± 0.9 

0.064 ± 0.001 



-504.2 ± 0.08 

-39.3 ± 0.1 

0.0799 ± 0.0002 



second layer 



Ea. (K) 

iE/N)o (K) 

{E/N)o - E^ (K) 

PO (A-^) 



-142.4 ± 0.5 

-166.8 ± 0.5 

-24.4 ± 0.7 

0.068 ± 0.002 



-178.9 ± 0.6 

-36.5 ± 0.8 

0.078 ± 0.003 



more so in the solid case, and the energy per D2 molecule 
differences between liquid and a solid phases are larger 
in the case of a second layer solid than for a single deu- 
terium sheet. The equilibrium density of the second layer 
solid is also comparable to the one extracted from the 
Maxwell construction for the entire system (0.078 versus 
0.070 A-2). 
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FIG. 3: (Color online) Same as in Fig. [T] but for D2 instead 
of H2. The lower limit for the surface density is fixed as the 
inverse of the experimental density for a promotion of D2 to 
a third layer (5.8 A^ Ref H). 
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FIG. 4: (Color online) Same as in Fig. [21 but for D2 instead 
ofH2. 



IV. CONCLUSIONS 

In this work, we have studied the complete phase dia- 
gram of the second layer of both H2 and D2 adsorbed on 
top of a single graphene layer. Since we used a set diffu- 
sion Monte Carlo calculations, the results correspond to 
the zero-temperature ground state of the system. To ob- 
tain the real stability limits, we had to perform Maxwell 
constructions between phases that comprised one and 
two hydrogen layers, but we found that if we used only 
the data corresponding to a second one, the description 
is good enough to reproduce the solid equilibrium of the 
complete two-sheet system. Our results for the promo- 
tion density to the second layer and the minimum den- 
sity in this second layer agree satisfactorily with available 
calorimetric results for graphite, pointing to the accu- 
racy of both the method used and the interaction po- 
tentials entering in the Hamiltonian. That our results 
on graphene are comparable to the experimental ones 
on graphite also means that both surfaces are basically 
equivalent as absorbents, being the only possible differ- 
ence the binding energy of the hydrogen molecules to the 
carbon surface. 

However, calorimetric measurements only give the to- 
tal density_for the onset of a solid structure at T — >■ 0. 
From Ref. 23, those densities appear to be 0.165 (H2) and 
0.178 (D2) A~^. Both are lower than the ones assigned in 
the same work to the oblique structures: 0.173 (H2) and 
0.186 (D2) A~^, respectively. The fact that the former 
densities are compatible with our simulation results sup- 
ports our suggestion of a double incommensurate solid 
as the structure for the inferred solid phase. Moreover, a 
comparison of the energies per molecule for the oblique 
solid and the two sets of triangular layers at the same den- 
sities indicates that, at least for T = 0, the commensu- 
rate arrangement is not stable. One could speculate that 
the disagreement between calorimetric data, in one side. 



and LEED and neutron diffraction, on the other, could 
be originated by the tiny difference between the diffrac- 
tion patterns of a double incommensurate structure and 
a double oblique layer i^^*^ In fact, the oblique phase was 
proposed in the past as the ground state of a second layer 
of *IIe on graphite in neutron diffraction studies^ but 
further calorimetric measurements concluded that trian- 
gular structures were preferred. This could be also the 
case for hydrogen. On the other hand, our results rely 
on empirical potentials that have been used to repro- 
duce reasonably well the equation of state of the first 
layer of both H2 and D2 on graphene;^ but we can not 
exclude that, in the future, more elaborate interactions 
could change our predictions for the second layers. 

On a different note, from the behavior of this second 
layer, we can say that when the number of adsorbed lay- 
ers grows, the entire arrangement becomes more "solid- 
like" than "liquid-like" . The reason is that, even though 
the binding energy in the second layer is lower than in 
the first one, the incommensurate solid structure is much 
more stabilized with respect to a liquid than in the first 
layer case. This is probably the reason why the experi- 
mental critical points for the liquid-vapor coexistence re- 
gions of II2 and D2 adsorbed on graphite approach those 
of the bulk solids as the number of adsorbed layers grows, 
^^ making impossible to obtain a stable liquid by reduc- 
ing the dimensionality of the system. 
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